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INTRODUCTION 

One  of  the  most  intriguing  operations  performed  in  various 
branches  of  Applied  Mathematics  is  the  transition  from  general 
equations  describing  the  state  of  a  spatially  distributed  material 
to  equations  applicable  to  materials  spread  out  in  a  relatively  thin 
layer.   An  example  from  Hydrodynamics  is  furnished  by  the  equations 
for  long  gravity  v/aves  in  shallow  v/ater  which  are  derived  from  the 
general  equations  for  irrotational  flow  of  an  incompressible 
liquid.   Usually  the  equations  resulting  from  such  a  metamorphosis 
bear  little  resemblance  to  the  original  ones.   Not  only  the  number 
of  the  independent  variables  but  also  the  order  and  type  may  be 
changed.   This  is  to  be  expected  from  the  fact  that  both  the  original 
equations  and  the  boundary  conditions  in  the  neglected  direction 
usually  contribute  the  new  equations.   The  present  report  attempts 
a  more  precise  analysis  of  the  transition  from  the  general 
3-dimensional  non-linear  equations  of  equilibrium  for  a  perfectly 

elastic  solid  to  the  equations  of  v.Karman  and  F8ppl  for  a  flat 

2 
plate  v;ith  force-free  faces. 

VJe  briefly  reviev;  the  situation.  We  deal  with  ^-dimensional 

transformations  x.  =  x.  +  u. ,  where  the  x.  are  the  cartesian 

ill  1 

coordinates  of  a  particle  in  the  natural  unstrained  state,  the  x. 
the  coordinates  of  the  same  particle  in  the  strained  state  and  the 
u.  the  displacements.   The  state  of  strain  is  described  by  the 
Quantities 


See  J.  J.  Stoker:  IJater  Waves,  Interscience  Publ.,  1957. 

o 

See  A.  F8ppl:   Vorlesungen  liber  theoretische  Mechanik,  v.  5  (1907); 
v.Karman,  Enc.  d.  math.  Wiss.,  v.  IV,  4  (1910);  W.  T.  Koiter:   Proc. 
of  the  lUTAM  Symposium  on  the  theory  of  thin  elastic  shells,  Delft, 
i960. 


employing  the  summation  convention  and  the  customary  notation  for 
partial  derivatives.   The  physical  characteristics  of  the  material 
are  given  by  the  strain  energy  density  function  \i ,    which  in  the 
case  of  an  isotropic  homogeneous  material  is  expressible  in  terms 

of  the  eigenvalue  of  the  matrix  of  the  e.,  .   For  small  strains  v;e 

^  xk 

have  up  to  terms  of  higher  order 

(^^)         ^'^  =  2(/+<r)  (^ik^ik  +  TTTo^^ii^kk^ 

(E  =  Young's  modulus,  (T=   Poisson's  ratio).   Suitable  stress 
components  t..  are  then  defined  by 

IT     \  4-        S'^ 

(Ic  )  t_. ,.  =  ^ 


ik   oe., 

ik 

The  equations  of  equilibrium  in  the  absence  of  body  forces  are 
then  given  by 


where 


(I«)  fjK-  2«0r*rk-^S_j^  . 


/-■ 


In  what  follows  we  assume  that  11  is  given  precisely  by  the 
expression  (lb)  (with  a  numerical  value  for  (T   between  0  and  -g )  > 
not  only  up  to  terms  of  higher  order.  This  special  choice  of  the 
strain  energy  function  is  not  at  all  essential  for  the  results  to 
follojw.   It  only  has  the  advantage  for  the  reader  to  give  the 
stress-strain  law  (Ic)  the  explicit  linear  form 


Our  results  vjill  not  depend  on  the  strict  linearity  of  these 

relations;  they  will,  on  the  other  hand,  only  be  valid  under  the 

assumption  that  the  stresses  or  strains  are  so  small  that  (If) 

represents  a  good  approximation  to  the  actual  stress-strain  lav;. 

The  P.,  then  become  quadratic  forms  in  the  t   .   We  shall  normalize 
ik  rs 

the  unit  of  force  in  such  a  way  that  the  Lame  shear  modulus 

-    E 
^   ~    2(1  +T) 

has  the  value  1.  This  has  the  effect  of  making  the  stress 
components  t..  dimensionless. 

Since  our  boundary  conditions  only  involve  stresses  it  is 
convenient  to  have  differential  equations  in  terms  of  stresses  only. 
We  then  have  to  supplement  the  equilibrium  equations  (Id)  by  the 
compatibility  conditions  which  express  the  vanishing  of  the  Riemann 
curvature  tensor  associated  with  the  transformation  x  =  x(x). 
Those  conditions  have  the  form 

^^S)   t^k,ij  "*■  ^ij,hk  -  ^hj,ik  -  ^ik.hj 


a 


(t_  ..5,,  +t_  K,5,  ,-t_  ,,5^,  -t_  ^,6,,  )  +q,,,  ,,,  =  0 


1+0"  ss,ij  hk   ss,hk  ij   ss,ik  hj   ss,hj  ik   ^hijk 

where  the  q,  .  .,  are  Quadratic  forms  in  the  first  derivatives  of  the 

^hi ,]  k 

t.,  with  coefficients  depending  on  the  t.,  .   Equations  (Id)  and 
(If)  form  an  overdetermined  system  of  mutually  dependent  differen- 
tial equations  for  the  symmetric  matrix  formed  by  the  t..  .   Under 
the  assumption  that  the  stresses  are  not  too  large,  say 


(Ih)  ^Ik^lk-^"' 

the  coefficients  of  the  highest  order  derivatives  in  the  equations 
will  not  differ  much  from  constants. 

In  the  case  of  a  "plate"  we  deal  with  a  cylindrical  region 
bounded  by  the  planes  x,  =  h  and  x^  =  -h  and  by  a  cylinder  with 
generators  parallel  to  the  x,-axis.  VJe  assume  that  no  forces  act 
on  the  plane  faces: 

(Ila)       t, -,  =  to^  =  t,^  =  0   for   X,  =  ±h  . 
Ip    25    35  3 

The  equations  of  v.Karman  and  F^ppl  represent  an  approximation 
which  connects  the  average  horizontal  stress  components  with  the 
deflection  of  the  middle  surface.   VJe  introduce  the  quantities 

h 

(lib)        \p(^l-^2)  -  h    f     ^*aP+Pap)^^3  '  ^""'^  ^   ^'^^ 

-h 

which  are  essentially  the  stress  components  t  g  averaged  over  the 
thickness.   On  the  basis  of  the  equilibrium  equations  (Id)  and 
boundary  conditions  (Ila)  the  to  can  be  expressed  with  the  help  of 
an  Airy  function  0  as 

(lie)        t^^  =  ^,22  '    h2  =  -^12  '        \2   =  ^,11  • 

We  also  introduce  the  quantity  w(x,,Xp)  =  u,(x-,Xp,0)  which 
represents  the  normal  displacement  of  the  middle  surface  x^  =  0. 
We  apply  a  rigid  motion  to  the  deformed  plate  to  give  it  approxi- 
mately the  same  position  as  the  original  one.   The  approximate 
equations  of  v.Karman  and  F8ppl  then  have  the  elegant  form 


(Ilia)  A^0  =  2(l  +  cr)(w^^2-^"^  11^  22^  ' 

(Illb)   A^w  =  i(.l-^)h-^(w^ll'0^22-2^^12^,12  +  '^22^,ll) 


where 


A  =  a^/ax^  +  h^/lx^   .      ^ 


Obviously  equations  (IIIa,b)  are  not  precise  consequences  of 
the  5-dimensional  equations  and  boundary  conditions.  Usually  one 
derives  them  on  the  basis  of  suitable  a  priori  assumptions  about 
the  relative  orders  of  magnitude  of  various  quantities.   Correction 
terms  can  be  obtained  by  looking  at  equations  (IIIa,b)  as  the 
lowest  order  result  in  an  asymptotic  expansion.  "^  Some  idea  about 
the  range  for  which  equations  (IIIa,b)  are  appropriate  can  be 
obtained  out  of  those  equations  themselves,  if  v;e  make  use  of  two 
"principles".  The  first  principle  asserts  that  in  an  approximate 
equation  all  terms  can  be  presumed  to  be  of  the  same  order;  other- 
wise the  equation  contains  superfluous  terms.   The  second  principle 
tells  us  that  the  order  of  magnitude  of  the  n-th  derivative  of  a 
solution  of  an  elliptic  system  of  partial  differential  equations 
can  be  obtained  by  dividing  the  order  of  the  solution  itself  by  a"^, 
where  a  is  the  distance  from  the  boundary  of  the  region  of  defini- 
tion (or,  for  smooth  boundary  conditions,  some  ether  typical  length 
associated  with  the  region).   Let  then  e  be  the  order  of  magnitude 


If  we  do  not  normalize  E/2(1  +  «j')  to  have  the  value  1,  v;e  have  to 
multiply  the  right-hand  side  of  (Ilia)  and  the  left-hand  side  of 
(Illb)  by  this  factor. 
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See  J.  L.  Synge  and  VJ.  Z.  Chien:  The  intrinsic  theory  of  elastic 

shells  and  plates.   Applied  Mathematics,  Th.  v.Karman  anniversary 

volume  (19^1),  pp.  103-120).  Also  P.  John:   On  finite  deformations 

of  an  elastic  isotropic  material.   New  York  University,  Institute 

of  Mathematical  Sciences,  Report  IMM-NYU  250  (1958). 


of  the  (dlmensionless)  stresses  and  let  us  consider  points  a 
distance  a  away  from  the  edges  of  the  plate.   The  second  deriva- 
tives of  the  Airy  function  0  as  average  stresses  can  also  be 
expected  to  be  of  order  s.  The  left-hand  side  of  (Ilia)  can  then 
be  expected  to  be  of  order  e/a  .  This  means  then  that  the  second 
derivatives  w"  of  w  should  be  of  order  /e/a  and  hence  A^w  of  order 
yi'/a'^.   Comparing  with  (Illb)  we  find  that  e  must  be  of  the  order 
h  /a  with  w"  of  order  h/a  ;  hence,  after  a  suitable  rigid  motion, 
w  and  h  would  be  of  the  same  order  in  a  region  of  size  a. 

Ideal,  from  the  point  of  view  of  applications,  v/ould  be  an 
estimate  for  the  difference  of  the  solution  of  equations  {Id,f ) 
with  boundary  conditions  (Ila)  supplemented  by  boundary  conditions 
on  the  lateral  surfaces  and  of  the  solution  of  equations  (IIIa,b) 

with  suitably  averaged  boundary  conditions  on  the  lateral  surface. 

2   2 
This  ideal  is  difficult  to  achieve.   IJhen  e  is  of  the  order  h  /a 

we  are  in  the  range  where  buckling  can  occur  and  the  dependence  of 

the  solution  on  lateral  data  becomes  irregular.   In  addition  the 

validity  of  the  v.Karman-F^ppl  equations  becomes  more  doubtful  as 

we  approach  the  edges.   Boundary  layers  might  develop.   It  seems 

safer  to  look  at  equations  (IIIa,b)  as  "interior  limits",  only 

valid  at  a  safe  distance  from  the  edges. 

In  the  present  investigation  we  consider  a  solution  t..  of 

the  exact  equations  (Id,f )  and  boundary  conditions  (Ila)  defined 

in  a  cylindrical  region  of  the  type  mentioned  and  satisfying  there 

the  inequality  (Ih)  with  a  certain  s.  The  actual  plate  itself  is 

thought  to  extend  laterally  beyond  the  cylinder  and  the  stress 

system  produced  in  the  cylinder  can  be  thought  of  as  produced  by 


appropriate  forces  applied  further  out,  \Je   shall  not  discuss  the 
delicate  question  of  what  kind  of  forces  actually  would  lead  to 
stresses  satisfying  (Ih)  in  the  cylinder.   VJe  shall  also  assume 
that  the  t.,  possess  derivatives  of  any  desired  order  in  the  closure 
of  the  cylinder,  but  shall  not  make  any  quantitative  assumptions 
whatsoever  on  those  derivatives.   Only  the  t.,  themselves  are 

restricted  quantitatively  by  the  a  priori  assumption  (Ih).   V/e 

2 
assume  that  e  lies  below  a  certain  universal  constant  e  to 

o 

guarantee  that  the  terms  v/ith  constant  coefficients  predominate  in 
the  equilibrium  equations  (Id).   VJe  consider  a  curcular  sub- 
cylinder  of  radius  2a  of  our  domain,  and  a  coaxial  cylinder  of 
radius  a,  where  a,  e,   h  shall  be  related  by 

(IVa)  e  <  h^/a^  <  1  . 

Then  in  that  smaller  sub-cylinder  the  v.Karman-Fflppl  equations 

2   ^' 
(IIIa,b)  are  valid  with  respective  maximum  errors  Keh  /a'  and 

Keh/a  ,  where  K  is  a  universal  numerical  constant.   N  tice  that 

the  errors  have  an  order  differing  from  the  order  that  is  to  be 

expected  for  the  terms  appearing  in  the  v.Karman-F^ppl  equations 

2   2 
by  the  factor  h  /a  ,  and  to  that  extent  can  be  considered  as 

negligible. 

The  error  estimate  derived  here  is  concrete  in  the  sense  that 

it  only  involves  a  universal  constant.   It  is  based  on  estimates 


It  is  likely  that  estimates  of  the  form  (Ih)  in  a  sub-cylinder 
can  be  derived  from  upper  bounds  for  the  total  potential  energy  of 
the  plate.   If  the  equilibrium  position  is  one  of  minimum  energy 
then  such  upper  bounds  can  be  obtained  from  trial  functions. 
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for  all  the  t..  and  their  derivatives  of  various  orders  that  are 
obtained  from  (Ih)  by  standard  arguments  used  in  the  theory  of 
elliptic  partial  differential  equations.   The  situation  is  somev;hat 
complicated  by  the  fact  that  boundary  conditions  for  only  three  of 
the  stress  components  t.,  are  prescribed  by  (Ila).   An  additional 
difficulty  for  thin  plates  arises  from  the  necessity  to  keep  the 
numbers  of  factors  h  in  the  denominators  to  a  minimum  if  adequate 
estimates  are  to  be  obtained.   It  is  to  be  expected  that  deriva- 
tives in  directions  parallel  to  the  faces  of  the  plate,  for  which, 
so  to  speak,  a  is  the  representative  length  will  be  of  smaller 
order  than  derivatives  normal  to  the  faces.   One  finds  indeed  in 
the  cylinder  of  radius  a  for  any  n-th  order  derivatives  of  a  t., 
involving  only  differentiations  with  respect  to  x,  and  Xp  the 
estimate 

(ivb)  \t[lh   1  V/^"" 

while  for  the  general  n-th  order  derivative 

(IVc)  \t[lU   1  V/a^-^h  . 

Here  C  is  a  numerical  constant  depending  only  on  n.  The  boundary 
conditions  (Ila)  then  imply  differences  in  orders  of  magnitudes  of 
different  stress  components  t.. ,  .   It  turns  out  that  In  the  cylinder 
of  radius  a 

(IVd)  It^^l  <  Deh^/a^  ,    It^^l  <  Deh/a 

with  a  certain  universal  constant  D.   This  represents  a  justifica- 
tion for  some  of  the  a  priori  assumptions  often  made  in  the 


derivation  of  equations  (IIIa,b).   On  the  basis  of  the  basis  of 
the  estimates  (IVb,c)  the  derivatives  of  the  normal  displacement 
w  can  then  be  estimated.   Here  a  preliminary  rigid  motion  has  to 
be  applied  to  the  deformed  plate,  to  make  the  transformation 
rotation-free  at  some  point  of  the  cylinder.   Once  all  quantities 
involved  in  the  3-diniensional  problem  are  estimated  properly  one 
obtains  immediately  the  previously  described  estimate  for  the 
errors  in  equations  (IIIa,b), 

The  main  labor  connected  with  establishing  this  result  lies 
in  the  proof  of  the  estimates  (IVb,c).   Here  the  principal  tools 
are  inequalities  between  Hilbert  norms  of  various  derivatives  and 
the  lemma  of  Sobolev  which  permits  transition  from  estimates  for 
the  norms  to  estimates  at  points.   There  is  nothing  new  in  these 
auxiliary  inequalities  except  the  special  care  that  had  to  be  taken 
in  bringing  out  the  precise  dependence  of  the  coefficients  they 
contain  on  a  and  h.   For  that  reason  short  proofs  of  the  general 
auxiliary  inequalities  have  been  included,  making  the  whole  report 
self-contained. 
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1.  Transformations 

VJe  refer  particles  in  their  unstrained  position  to  cartesian 
coordinates  x, ,  x^,  x^.   The  same  particles  in  their  deformed 
position  shall  have  cartesian  coordinates  x.  =  x. (x, ,Xp,x_) . 
Generally  for  a  function  f{x.,x^,x^)   vie   denote  by  f  ..  the  partial 
derivative  3f/cx. .   We  put 

(1-1)  Pik  =^i,k  • 

The  p.-.    are  then  the  elements  of  "Jacobian  matrix".   They  satisfy 

We  assume  for  the  Jacobian  determinant 


do)  det  (p^j^)  >  0  . 

We  can  then  introduce  the  p"""  that  form  the  reciprocal  of  the 
Jacobian  matrix  by 

(1.4)  p'Sk  =  PljP'""  '  ^Ik 

(using  the  summation  convention). 

The  line  element  ds"  in  the  deformed  state  is  given  by 

(1.5)  ds"  =  dx\dx^  =  °iic^''^i^^k 


whe 


re  the  g..  are  the  components  of  the  metric  tensor  given  by 


(1.6)  '  Sik  =  PjiPjk  =  Ski  • 


The  reciprocal  matrix  is  defined  by 


11 


(1.7)  g^^gjj,  =  S.jgJ^  =  5,^ 
and  also  given  by 

(1.8)  gi^  =  piJp^^'  =  s^i  . 
Customarily  the  quantities 


(1.9)  ^ik  =  |(gik-^ik) 


are  considered  as  the  components  of  the  strain  tensor. 

We  shall  make  use  of  the  Christoffel  symbols  of  the  first  kind 

(1.10)  rK,iji  =i(gj^_j^H-g,^^_j-g.j_^) 

and  of  the  second  kind 

(1.11)  [i^J  =  g^"[r,i,]  . 

They  are  connected  with  the  derivatives  of  the  Jacobian  matrix  by 
the  relations 

(1.12)  Pgj^^  =  P^^[k,ij] 

(1.13)  P'^,  =  -P^"[r^\- 

As  a  consequence  of  the  fact  that  we  are  dealing  with  trans- 
formations from  euclidean  3-space  into  itself  the  g.,  satisfy  the 
identities  arising  from  the  vanishing  of  the  Riemann  curvature 
tensor: 

°  "  ^ijk  "  2  ^%k,ij  "^^ijjhk  "Shj,ik' Sik,hj^ 


-.  [a,hk]  [^.M  -  [a,hj]  i.\]  . 
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2.  Stresses  and  equilibrium 

The  potential  energy  In  a  region  is  given  by 


/ 


W  dx 


where  dx  =  dx,dx^dx^.   VJe  assume  that  the  medium  is  homogeneous  and 
isotropic.   For  simplicity  we  take  the  case  where  the  s train-energy - 
density  W  is  given  by  the  standard  expression 

(2.1)  ^-^  =  t^^ik^k  -^  I  ^ii^kk 

where  the  e..  are  defined  by  (1.9)  (in  contrast  to  other  definitions 
for  the  e.,  that  have  been  used  in  other  connections)  and  A,  |i  are 
the  Lame  constants.   Instead  of  A  and  p,  we  can  use  Young '  s  modulus 
E  and  Poisson's  ratio  or  where 

^^•^^        ^  =  (l+a-)(l-27)  '   ^   =  2(1+0-)   • 


We  assume  that 

11  >  n     n   <    -Y  <  ^ 

2 


(2.3)  p.  >  0  ,   0  <  r<  i  . 


We  can  express  W  from  (2.1),  (1.9)^  (1.6)  in  terms  of  the  p^^^ 
and  Introduce  the  quantities 

(2.^)       q.,  =  ^-^^  =  2LLp.  (e  ,  +  T-^-n^   e.  .5  ,  )  , 
^    '       ^ik   cip.,     '^^ir  rk  1  -  2<y^  jj  rk 

The  equations  of  equilibrium  in  the  absence  of  body  forces  then 
have  the  form 

(2.5)  Qik^k  =  °  • 
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The  force  acting  on  a  surface  element  in  the  deformed  position  has 
components 

(2.6)  q^^£,^   dS 

where  dS  is  the  area  of  the  surface  element  in  the  unstrained 
position  and  the  ^,  the  direction  cosines  of  the  normal  of  the 
surface  element  in  the  unstrained  position. 

VJe  introduce  a  dimensionless  strain  energy  density 

(2.7)  V  =  iw=e.,e.,  +  ,  ^A^  e.  .e,  , 

'  IJ,      ik  ik   1  -  ^T  11  kk 

and  corresponding  dimensionless  stress  components 

SV    1  cW 


(2.8)  t 


ik   oe . ,    u,  de , , 
xk   ^   ik 


(We  could  instead  choose  the  units  of  force  in  such  a  way  that 
M,  =  1.)   In  contrast  to  the  matrix  q.,  the  m.atrix  of  the  t.,  is 
symmetric: 

(2.9)  t,j^  =  t^i  . 

By  (2.8),  (1.9)  the  t.,  and  g.,  are  connected  by  the  formulae 

(2.10)  Sik-^k=  ^k  -TTF^jO^k 
while  by  (2.4),  (2.8) 

(2.11)  -.    q^i^  =  ^^Plr^k  • 


14 


The  equations  of  equilibrium  (2.5)    can  be  v;ritten  in  terms 
of  the   stresses   t^^.      From  (2.5),    (1.6),    (2.7) 

°  =  ^  Pij^ik,k  =  il  ^Pij^ik^k  -  il  Pij,k^ik 

_   /  a.     X  1  SW 

-  ^PijPir^rk\k  ~  il  Pik,j  "Sp" 


Ik 


=  (Sjrt^l,)^j,  -  V^j 


t  .,  ,  +  (2e  .  t  ,  )  ,  -  V 
jk,k      jr  rk',k    ,j 


or 


(2.12)  0  =  t .,  ,  +P  .,  , 

jk,k   jk,k 

where 

(2.13)  P  ,  =  2e  .  t  ,  -  V5., 

jk     jr  rk     jk 

"  ^jr^rk  "  T+T  ^ii^Jk 

"  ¥  ^^ir^ir  "  r+7  ^ii^rr^^jk  * 

Equations  (2.12)  represent  three  first  order  partial  differential 
equations  for  the  quantities  t  ,  .   In  addition,  the  t.,  satisfy  the 

J  K  IK 

compatibility  conditions  arising  by  writing  equations  (l.l4)  in 

terms  of  the  t.,  instead  of  the  g.,  .   Putting  T  =  t.  .  we  have 
ik  ^ik        -^      JO 

(2.14)  0  =  t.  ,  ,  .  +  t.  .  .  ,  -  t,  .  .,  -  t.,  ,  . 

hk,ij   ij,hk   hj,ik   ik,hj 

-  TT^  ^^i/hk-^  ^hk^io  -  ^,ik^hj  -  ^,hj^ik)  -^^hijk 
where 
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There  are  essentially  six  compatibility  conditions.   VJritten 
out  in  detail  they  look  like  this: 

(2.l6a)    ^11^22  "^^22, 11  ~  ^^12,12  '  TTa"  ^^,11  "^^,22^  "  ^1212 

(2.l6b)         t22^33+t55^22  "2^23,23   "   1+7   ^^^,22  "^  ^,  33^    "    ^2323 
(2.16c)         t33^,^+t^3^^-^-2t^^^^-^    __J_^(T^^^+T^^3_)    =   q.^3^ 

(2.l6d)         ^23, 12  ^^12,23 '^22,13"  ^13,22   "^  m?^, 13     "^2132 
(2.l6e)       t^^^ 23 +^23, 31 '^33, 21 -■'^21,33  "^  rr:? '^,21     "^3213 


(2.l6f)       ^12,31  "^^31,12  "■'^11,32  '^32,11  "*"  TT^  ^,32     "^1321 


Putting  o=i,  k=hin(2.l4)  and  using  (2.12)  we  get 


(2.17a)   0  =  T^^^  +  ^   (P^,^,^  +  I  q^iij,)  . 


Just  putting  i  =  j  in  (2.14)  then  yields 


(2.17b)    0  =  t^j^^,,  +  j^   T^^^  +  ^  (P-i^,.  +  2  ^jiij^^k 

"^  ^hi ,  ik  "^  \l ,  ih '*"*^hiik 
or 

where 


16 


<2-W)   S^^.  -T^ffjl.lJ-^I^Jllj'Shk-fhl.lk-^l.lh-^ilk 

For  k  =  h  we  have  then 


(2.17e)   T  ..  =  4^S,. 
f  11    d  +  ry     11 


5.  General  Inequalities  needed 

VJe  consider  tensors  u  with  components  indicated  by  subscripts. 
For  example,  u  =  ^^iikaS^*  ^^^^   Roman  subscripts  run  through  1,  2, 
5,  Greek  subscripts  take  the  values  1  and  2.   A  comma  preceding  a 
subscript  1  or  a  indicates  partial  differentiation  with  respect  to 
the  corresponding  variable  x. ,  respectively  x  .   No  laws  of  trans- 
formation  for  tensors  will  be  invoked. 

For  tensors  we  have  the  operation  of  multiplication.   If,  for 
example,  u  =  (u.),  v  =  ("^ij^^  then  uv  -  w  =  (w^^^)  where 

w.  .,  =  U.V.,  .   Addition  of  tensors  of  the  same  kind  and  multiplica- 
ijk    1  jk 

tion  by  a  scalar  are  defined  as  usual.  Accordingly  the  square  of  a 
tensor  u  =  (u.,  )  is  defined  by 

u^  =  uu  =  (u,.j^u.^)  . 
The  length  |u|  of  a  tensor  u  =  (u.  .)  shall  be  given  by 

(3.1)  |u|2  =  u.^u.j^  , 

using  the  summation  convention.   Then 

(3.2)  |uv|  =  lu! |v|  ,    |u+v|  <  |u|  +  |vl  . 
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Differentiation  of  tensors  leads  to  new  tensors.  For  u  =  (u,,  )  we 
write  u'  =  (u^^  .),  u  =  (u.,   ).  Then  the  usual  rules  for 
differentiation  of  products  apply. 

We  shall  consider  u  defined  in  a  slab  -h  j^  x^  _^  h.   We  write 


/  udx   for    ///  udx-dXpdx 


|Xj|<h 


2--3  • 
The  norm  of  a  tensor  u  is  defined  by 

(3.3)  Hull  =/  f  |u|2  dx  . 

For  quantities  X,  Y  we  use  the  symbolic  relation 

X  =  0(Y) 

to  indicate  that  there  exists  a  universal  constant  K  with  a  positive 
numerical  value  such  that 

|X|  <  KY  , 

where  K  is  computable  in  principle  in  terms  of  universal  constants 
introduced  previously.   Similarly  X  =  o(Y)  means  that  |X|  ^  KY 
where  K  can  be  made  smaller  than  any  preassigned  universal  constant. 
Of  course  special  care  must  be  taken  in  concluding  that 
0(o(Y))  =  o{Y);  we  shall  only  make  use  of  this  relation  where  it  is 
clear  from  the  context  that  the  constants  indicated  by  "0"  and  "o" 
do  not  depend  on  each  other.  We  shall  often  have  occasion  to  use 
the  estimate 

(3.^)  luvl  =  O(lul^)  +  o(|v|^) 

which  follows  from  the  inequality 
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2luv|  <  C|u|^+C-^|v|2 


which  is  valid  for  any  positive  C. 

We  Introduce  an  auxiliary  function  C  depending  on  a  point 
(x,,Xp)  and  on  a  positive  parameter  a  by 

/^(l-a"^r^)^   for   r  <  a 

(3.5a)  C  =  ) 

(  0  for   r  >  a 

where 


(3.5b)       r  =J(x3^-x°)2  +  (x2-x°)2  . 

The  function  ^  is  of  compact  support  in  x, ,  Xp,  non-negative,  and 
has  continuous  first  derivatives.  Moreover 

(3.5c)       kl  =  0(1)  ,   \cU  0(a-^)  ,   ICi  =  0(a-^)  , 

(3.5d)      \t\^   =  O(crCl)  . 

For  any  scalar  u(x^,X2,x,) 

(3.6a)  He  u||  =  0(Hl'  u||)   +  o(|tC  u  ^^11) 

(3.6b)  lU  u|l=  0(|K  ^^aJ-'  11*^*^1'^ 

(3.6c)  IK  u^ll  =  OdCuail  +  llC  u^ll) 

(3,6d)       IK  u||  =  0(1|C  utl  +  A-^IU  u||)  +  o(Al|Cu||)  . 
Proof:   By  integration  by  parts,  using  (3.5d), 
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=  o(f  {\'i  u|K  u|+  Ul^luilu^ppDdxJ 

=  0(||C  u||2)+o(|K  u        ||2+  11^  u||2) 

PP 

from  which  (3.6a)  follo\\rs.   Similarly  (3.6b)  follows  from  the 
identity 

lU  'ull^  =jr^2  u^^pu^^p  dx 

=  0(|U  u^aall^+li*^*  u|h+o(|K  uf)    , 
(3.6c)   from  the  identity 

IICu^f.y^rA^u^^u^pU^pdx 


o(/(lCuu|K  u2|+CKI|up|u|)dxJ) 


=  0(|kuul|2+  lie  u||2)  +o(||C  n^f)    , 


(3.6d)   from 
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Cut 


^  ^,a  ^,a  ^^ 


More  complicated  to  derive  are  estimates  involving  derivatives 
with  respect  to  x,.  Let  f(s)  denote  a  complex  valued  function 
defined  for  |s|  ^  h.  Then  for  any  real  i 

h  .    h 

(3.7)   /  ^^If'l^ds  =  oy  ((4Vh-2i^2j|^|2^|^„|2j^\ 

-h  -h  ^ 

Proof:  Let  first  -h  <  s  <  0  and  0  <  !;"■'•  <  h,  so  that 


Then 


-1^' 


lr-'f'(s)| 


-h  <  s+^"-^  <  h   . 


s-.r^ 


f(s+r-^) -f(s)  -J    (s+r^-t)f"{t) 


dt 


s+i 


-1 


=  o(|f(s+r^)i^+  |f(s)|2+r^/    if"(t)i^dt) 


s 


Integrating  with  respect  to  s  from  -h  to  0  yields 


h 


/  r^if'i 


r^if'l^ds  =  0   /    (if|2  +  r'^if"i^)ds'    . 


-h 


-h 


y 


-1  2 

The  same  inequality  holds  for  -h  ^  i~     <  0,    since  only  i     enters, 

Replacing  s  by  -s  and  adding  the  resulting  inequalities  we  find 
that  for  0  <  Ul"-"-  <  h 


/   , 


(..:)-i! 
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h  , 

i 

-h  -h 


r^lf'l^ds  =   o(    f     (|f|2+r^|f"|^)ds\ 


^■1-1 


If    Ur     >  h  we  find  then 


h 


J   r^if'i^ds  =  rV^ j   h^if'i^ds 

-h  ^-h 


=  0(  rV^  f     {|f|2  +  h^f"|^)ds) 

"^       -h  ^ 

/  h  ^  . 

=  on    (e-V^ifi^+rNflMsi  . 


-h 


For  either  range  of  |sl~  the  inequality  (3.7)  follows  by 

multiplication  by  t,   , 

Replacing  i   by  ^/A  in  (3.7)  we  get  the  formally  more  general 

inequality 

h  /-  h 

(3.8)  f     i^\f'\^dB   =   0{^f     ((^^A-2  +  ^V2)|f|2  +  A2|f"|2)ds\ 

-h  -h  y 

valid  for  complex  values  f(s)  and  real  i,   A.  For  any  scalar 
u(x,,Xp,x,)  we  can  apply  (3.S)  to  the  function  f(s)  =  u(x-,X2,s). 

Taking  here  4  =  C~  IClh"  ,  multiplying  the  resulting  inequality  by 

2  2 
C  h  and  integrating  with  respect  to  x-^,   Xg,  we  obtain,  using 

(3.5d), 

(3.9)  Ik  u^^ll  =  0(h-^A-liU*u!|+h-^iKu|l+Ah|Uu^^^||)  . 

Let  u(^,Xp,x,)  denote  the  Fourier  transform  of  u(x,,X2,x,) 
with  respect  to  x,.   VJe  apply  (3.S)  to  f(s)  =  u(^,X2,s). 
Integrating  with  respect  to  4  and  x^   we  find  from  Parseval's 
equation 


22 


J   |u  ^^I'dx  =   /    ii^\u  ,1^  d^dx^dx 


2,  .         /    ,.2|A      ,2 

,3'      -— 2^'^5 


=   o(j    ((e'^A-2  +  eV2)|i|2+A2lu^^3l2)   dldxgdx^^ 

A   similar  estimate  holds  for  u   ^r>-      It  follows  that 

\\nj\  =   0(A-ll|ull  +h-l|lu||+A||u^33ll)    . 

Replacing  here  u  by   Cu,    using   (^•9)    with  A   replaced  by  A/h  and 
(3»6a)   we  find 

(3.10)        lU   u^jll  =   0(A-1(||C   u||  +  Hrul|)+A|U  u^^^ll 


+ 


h-i(||c  ull  +  iu  uID). 


This  formula  gives  an  estimate  for  the  mixed  derivatives  u  ,  in 


terms  of  "pure"  derivatives  of  u. 


Let  again  f(s)  be  defined  for  |s|  j^  h  and  let  F  =  Max  |f(s) 

|s|<h 
Then 

h  /  h  . 

(3.11)   f    r'^{s)ds  =  o{^J    (p2f"2(s) +h-V|dsj. 

'-h  -h  ^ 

Proof:  We  first  consider  the  case  where  f(±h)  =  0.  Then 

h  h 

J    f'^  ds  =  -4  J  f'2f>'f  ds 
-h  -h 

=  off    f2f"2  dsj  -^  oQ'    f'^  ds) 

-h        -^      -h      ^ 
and  hence 
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P      ^'^  A„    -    nf    P     ^2,"2  .J^ 


f'"  ds  =  0(  /   f^f"^  ds'   . 


-h  ^  -h 

Let  next  f(s)  be  arbitrary.   We  introduce  the  linear  function 

g(s)  =  f'"'^-/'-^'  s  +  f'h)^/'-h) 

SO   that  f-g  =   0  for  s  =   ±h.      If  F  =  Maximum   jfj    we  have    |g|    =   0(F) 
and    |g'  I    =   OlFh"-"-).      Then 


J    f'^as  =  o(^f    ((f'-g-)^  +g'^)ds^ 


-h  '  -h 

h 


=  oQ     ((f"-g")^(f-g)2+g''^)ds) 


•h 
=   off     (f"^(f^+g")+g'")dsj 


."2,^2,_2,    ,      A       ^ 


=   o(  f     (p2f"2+h-V)ds)      . 


Taking  for  f(s)  in  (3'11)  a  function  u(x,,Xp,s)  where  |ui  _^  U 

2 
in  the  support  of  C  we  find  multiplying  by  C  and  integrating  with 


respect  to  x,  and  Xp  that 


(3.12)       lU  u^^ll  =  0(Ul|Cu^^^||  +  |Uh-V||)  . 


-  c- J 
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Sobolev's  Inequality  for  a  slab 

At  any  point  x  =  (x,  ,x^.,x^)  of  the  slab  Ix^l  <  h 

jt   d   ^  J?  — 

(5.13)   h^/2  ^-1|^  ^1  =o(N(llCu||  +  irC  ulD  +  a-^iKull 


+  N-l(h2a-2||Cu^3^||+a-2||Cu||)) 


where 


N  =  jTThV^  . 


Proof: 


For  X  outside  the  support  of  ^   (3.13)  is  satisfied  trivially. 
Let  then  x  lie  in  the  support  of  C,   and  assume  without  restriction 
of  generality  that  -h  ^  x,  j^  0.   Consider  the  points  y  on  a  ray 
from  x: 

y  =  X  +s|  ,    s  ^  0  , 

the  direction  of  the  ray  being  determined  by  the  unit  vector 
i  =  (|,,^p,^,).  VJe  restrict  the  direction  of  the  ray  by  the 
assumption  that 

(3.14)  0-  S'E  ■ 

For  u  defined  in  the  slab  the  expression  C(y)u{y)  =  ^(x+s^)u(x+sC) 
as  a  function  of  s  will  then  be  defined  for  all  s  ^  0  and  have  its 
support  in  the  interval  0  j^  s  _^  Na.   Indeed  we  have  for  0  _^  s  _^  Na 

-h  <  x^  <  y^  =  x^+s'^^  1  sl^  <  h 


See  S.  Sobolev,  Mat.  Sbornik,  N.S.  4  (1938),  pp.  471-497; 
K.  0.  Friedrichs,  Comm.  Pure  Appl.  Math.  6(1953),  PP.  29^-325; 
L.  Nirenberg,  Comm.  Pure  Appl.  Math.  9  (1956),  pp.  50^^-530. 


25 
so  that  y  lies   in  the   slab.      For   s  =*  Na 

7(y^-x^)'^  +  {y2-x^)''^  >  J  (yi-x^)''^+(y2-X2)^   -  J{x^-x°)'^+(x2-x°)^ 

>   s  /if"+~i|-  a  =   s  Jl-  d  -a 


>  JN^aTh^  -a  =  a 


so  that  C(y)u{y)   =  0.      Then 


^    00         2 
^^(x)u^(x)   =(  j      s  ^  C(x+s|)u(x+sl)ds' 


/CO  .2 


0 

00 


l^af    s^((C(y)u(y))^^j^  4^ej^)^ds 


Since   the   i     are  bounded  and   i-,   satisfies    (3.14)    we  have   then 
a  p 

/-co 

C^(x)u2(x)   =   0  Ua  j      (s^dC  u|2+  IC  u|+  IC  u|^) 

0 

+    sVN-V^dC  u^^l^+U  u^3|2)  +  s2hVV^|C  ^jjI'HsV 

Integrating  over  all  unit  vectors   i   satisfying   (3.14)    it  follows 
that 


•  9  •   • 


hl/2  a-l|C  u|    =  0(N(|Ku|l  +  |U  ui+llC  u||) 


+  ha-l(|lC  u^^lU   IK  u^^lD+hW^IICu^^^ll)    . 
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-l^T-1 


If  we  use  here  (5.6a),  (3.9)  with  A  =  a'-'-N"  ,  (3.10)  with 

A  =  ha"  N'  and  (3.6d)  with  A  =  Na  the  inequality  takes  the  form 

(3.13). 

All  the  inequalities  derived  here  for  a  scalar  u  can  be 
extended  immediately  to  tensors  u  with  a  fixed  number  of  components. 

4.  Estimates  for  the  stresses  and  their  derivatives 

Assume  that  the  t,,  are  defined  and  satisfy  the  differential 
equations  (2.12),  (2.1^^)  for 

(4.1)  Ix^l  Ih  ,    (x3^-x°)2  +  (x2-x°)^  <  a^  . 
Assume  also  that  the  t-,,  satisfy  the  boundary  conditions 

3K 

(4.2)  t^^   =0   for   x^  =  ±h  ,        k  =  1,2,3 
and  that  the  tensor  t  with  components  t.,  satisfies 

(4.3)  |t|  <  6 

at  every  point  of  (4.1).   The  cylindrical  region  (4.1)  shall  have 
proportions  related  to  the  upper  bound  e  for  the  stresses  by  the 
inequalities 

(^•^)  ^i%4'  111- 

a 

Here  e  denotes  a  universal  constant  which  will  be  chosen 
o 

appropriately  small  in  the  arguments  to  follow.   In  order  that 
there  are  quantities  a,  h,  e  satisfying  (4.4)  at  all  we  must,  of 
course,  have 

(4.4a)  ^  ^  ^o  • 
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If  the  second  inequality  in  (4.4)  but  not  the  first  one  is  valid, 
we  can  always  cut  down  the  radius  a  of  the  cylinder  sufficiently 
to  satisfy  both  inequalities,  provided  (4.4a)  holds. 
Defining  C  as  in  (3.5a)  we  put 

(4.6a)  B  =  lUt"ll  , 

(4.6b)  A  =  Max  (lUt^p^^.ll,  lUt^p^^jll,  IUt,:5,piL  lUtj^^iJD 

ct,  p ,  / 

5,i,k 
so  that 

(4.6c)  A  j^  B  .' 

By  (2.12) 

IUt^3,5,ll  =  o(|Kt3„^^3linUp"ll)  , 

by  (2.17b)  with  h  =  k  =  3 

IICT^^^I!  =  odUt^^^^^ll  +  IUt^^^^^ll  +IUp"II+  llCqIl  )  , 

by  (2.17a),    (3.6b) 

H  n  -  0(|U  T^^^ll+IIC  Til) 

=  o(jU  T^^^ll  +|Up"II+IK  qll+HV  til)  , 

by  (2,l6b,c,e) 

+  IUt3„,p3ll+IICt5p_„3ll  +  lUql|)   , 
by   (2.17b)   with  h  =  a,   k  =  p  and   (3.6b) 
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Altogether  then 

To  estimate  the  right-hand  side  we  use  the  equations  {2.17c)  and 
(2.12),  observing  that  by  definition  |S|  =  0(|p"|  +|q|).  We  have 
by  integration  by  parts,  on  the  basis  of  the  boundary  conditions 
(4.2), 

ll^*3k,iall  ll«*3k,lall  '  /  «Sk,la*5k,ic  <»^ 

=  -/  '?'t5k,lc.lSk,a+2«,pbk,Pat5k,a'^== 
=  /('3k,ii(«Sk,a',a-2«,pV,aa*5k,a''S^ 
=  /'-TT^^3k'«Sk,a',.-V«Sk.al, 


a 


-  2«,p*5k,PaSk,a'«^ 


=  /'-rb^,5a'?Sk,a',k+V^Sk,a', 


a 


-  2??,e*3k,Pc.*5k,a''i^ 


=  /  'tT7  «'^,3a^3k,ka  +  IT?  2«,3^at33,a5^  ':2Skbk,aa 


+  2«,„S3^t5,^_„  -2«^pt3^^g„tj^^„)dx 
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=   0(B||C   P'il  +  liCP"f  +    IKq||^+    Ik   tf) 

=   0(aV2llCp'||2+    IkP-ll^  +  IICqll^  +  HC   tf)+o(hW+A2) 

where  we  have   made  use  of    (5.6a)    to  estimate    ||C  t||.      It  follows 
that 

(4.7)        A  =  0(ah-l||C  p'||  +  IUp"II+  lUqII  +  \U  1 1| )  +  o(ha-^B)    . 

By  definition  of  B 

B  =  0(A+    Max  (|kt„p_,^IUIUt,^_53ll)). 

Here  by  (2.17b)   with  h  =   3,   k  =  7 

Ut^y^^^\\=   0(A+   lUt„^^3y||+IKP"ll+IICqll)   . 

Using   (3.10)   with  A  replaced  by  1  and   (5»6d)   with  A  replaced  by  _a^ 
we  find  that 

=  0(A  +  ||V  til  +  h"^|U  tll  +  a'^h"^||C  t|l  +  ah"^A)   . 
Since  1  <  a/h  we  have  then 

B  =  0(ah-lA+irC   t|l+h-^iK  tll+a-V^llC  t||  +  |Up"||+  |Kql|  )  . 

Comparison  with   (4.7)    then  yields 
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(4.8a)  A  =   0(ah-lllCP'||+   ||CP"il  +  lUqII  +  11?  t|| 

+  a-llKtll+a-^llCtll) 

{4.8b)  B  =  0(aV2|Up'||+ah-^IU  P"ll  +  ah-^|Uqll 

+  ah-l|rCt||  +h-l|Kt||+  a-VllUt||)  . 

So  far  no  use  has  been  made  of  the  special  form  of  the 
expressions  P  and  q  and  of  the  inequality  (4.3).   We  observe  that 
by  (2.13)  the  components  of  P  are  quadratic  forms  in  the  t.,  with 
numerical  coefficients.  By  (2.15)  the  components  of  q  are  quadratic 
forms  in  the  first  derivatives  of  the  t. .  with  coefficients  that 
are  rational  in  the  t.,  and  bounded  for  |t|  ^  e  _^  e^.   It  follows 
that 

|q|  =  O(lfl^),   |P"I  =  O'itllfl  +lt'l2)  =  0(e|t"|+  tt'l^)  , 

IP'I    =  0(|t||t'l+  Itiit'l)   =   0(elt"|+   ItUt'l)    . 
Thus 

IUP"ll  +  IICqll  =   0(£|Kt"||+IUt'2i|)    =  0(eB+|Ut'2||) 
lie  P'll  =  0(eB+  lU   i2||^/^lkt'2||^/2j    ^ 

Here  by   (3.6c),    (3-12) 

(4.oa)       lU   t^ll  =  OdlCt  til  +IIC  t2|t)=0(eA+eV2|UII) 


(4.9b)    iut'^ii  =  o(tu  i^ii  +  iict^^ir ) 


/'    ^^^--Idi.^d 


2.-2. 


=  0(eA+eV'=^||Cll+  eJUt  „||  +  eV'^||Cll) 


=  0(eB+e''h-^||Cll)    . 


.    t 


31 

Then 

lU  P'll     =   0(sB  +  ah-^|K  t^llH-   ha-^IUt'^ll) 

=  0(eB  +  eah"-'-A+ e^a'-'-h'-'-llCll  )    • 
If  we  observe   that  also 

lU  til  =  0(e|Ul|)   =  0(£a-^llCll) 
JlY  tll=  OCelJCll)   =  0(£a-2|lcl|) 
we  find  from  (4.8a,b)   that 

A  =  0((l  +  ah"-^)£B+ea^h"^A+ (e^h"^+ ea"^)|KII  ) 
B  =   0((a^h"^+ah"-'-)eB +ea\"\+  (ea'^h'-^-H-  e^ah"^)||Cll  )    . 
Using   (4.^)   with  e     sufficiently  small  we  find  then  that 

A  =   0(ah"-^eB  +  ea"^IUll  )    ,        B  =  0(  ea^h'^A  +  ea"-'-h"-'-|U  II  ) 
and  hence  that 

(4.10)  B  =   0(£a"^h"^lU||  )    ,  A  =  OCea'^lUII)    • 
It  follows  from  (4.9a,b)   that 

(4.11)  lU  t2|l=  ©(eV^llCll)    ,  lUt"2|l  =  0(eV2|Ull)    . 
Recalling  the  definitions  of  A,    B  we   i«?ite   (4.10)    in  the   form 
(4.12a)       lUf'll  =   O(ea-V^lUII) 

(4.12b)       IK  tlMICt^p^^^ll,  lU  t^^lMlCt^^ll  =  0(ea-2lK||)    . 
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For  the  purpose  of  deriving  analogous  inequalities  for  the 
norms  of  the  higher  order  derivatives  of  t  we  have  to  vary  the 
parameter  a  entering  the  definition  of  the  auxiliary  function  Cj  we 
write  then  everyv;here  C,^   for  the  function  denoted  up  to  now  by  C- 
The  parameters  x^  and  y^  will  be  held  fixed.   We  write  then  (4.8b) 
in  the  form 


(4.13a)    Ik^^f'll  -  0(aV^  lU^P  '  11+ ah-^HC^P"  li  +  ah-^IU^  ql| 

+  ah-^  ir4t||+  h-1  lU^tll  +  a-VMU^tll 

and  a  special  case  of  (4.3a)  as 

(4.13b)    lU^tll  =  0(ah-l  lU^r'H  +  IIC^P-ll+llC^qll 

4k^t\\+a-^   IUat||  +  a-2  lU^tll  )  . 

In  deriving  these  inequalities  we  only  made  use  of  the  differential 
equations  (2.12),  (2.14),  the  boundary  conditions  (4.2),  and  finally 
of  the  fact  that  ha~  ^   1,  or,  rather  only  that  ha"  =  0(1).   The 
special  form  of  the  expressions  P  and  q  and  the  limitation  (4.3)  on 
t  were  not  used.   It  is  clear  then  that  we  can  replace  everywhere 
in  (4.l5a,b)  t  by  t,  if  we  replace  at  the  same  time  P  and  q  by  P 
and  q,  since  t  satisfies  the  same  kind  of  differential  equations 
and  boundary  conditions.   We  can  at  the  same  time  replace  a  by  a/2. 
The  resulting  formulae  become  simpler  if  we  observe  that 

Finally,  we  make  use  of  the  fact  that  as  a  consequence  of  the 
differential  equations  (2.12),  (2.14) 
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|t'"|  =  o(|t  |+|p'"l+  k'l)  . 

We  then  obtain  the  inequalities 

(4.14b)    llC^/gf"!!  =  0(aV2  IK^/gV'll+ah-^  Ue,/2^'''\\ 

H-ah-^  lk^/2^;!l+a-Vl  IIC^^II) 

(4.14b)    llC^^/gt'll  =  0(ah-l  IUa/2'^'11  +  II  V2^'"ll 

+  IUa/2q'll  +  a-2|U^i||)   f 
Put 


Then  ^ 

lP"'|  +  |q'l    =   0(|tt'"l+  |t't"|+  It'^!) 

=  o(e|t'"|  +  |t'||t""|  +  |t'Mf2!) 

Ip'l    =  0(|t  t"|+  It  t"|  +  It'  *t|)   =  o(£|t"'|  +  |t||t"|+  It'Mtj)    . 
It  follows  then  from  (4.11),    (4.12a,b)  using  ^    /g  =  ^(^a/2^a^   ^^^* 
11^2^ "'II -^11^2^ 'II  =  °(^  IK^/2t"^||+M^ea-V^  lU^JI 

|U^/2^'ll  =   0(e  ||C^/2^'"II  +  em3_a-Vl  ||C^|| +M^ea-^||)    . 

Substituting  in  (4.l4a,b)  and  making  use  of  (4.4)  results  for 
sufficiently  small  e^  in  the  relations 

11^2^'""  =  0(eah-\  ||Ca||  +  £h-\  Ik^ll+a-V^  IK^tJI) 


i^ii^;!  '  ;■' 


;V  *j|  +  . 
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+  ea-V\  lK^II+a-2  ||^^t||) 

=  0(eh-2m^  IUall+£a-V\  11^^11+ a'^  JU^ill  )  . 
Here  by   (4.11) 

(4.15a)  lU^tll  =   0(|U^l|l/2  ||C^;2,|l/2j   ^  Q(^^-l  11^^,1) 

(l|.15b)  llC^t'll  =   0(|UJ|^/2  lU^f  2l|l/2)   ^   Q(^j^-1  11^^,1)    ^ 

VJe  now  make  use  of  Sobolev's  lemma   (3. 13).     For  us 

2   -2 
2  _^  N  j^  5  since  ha       <  1.     We  observe  that 

IIC^ll  =  0(ah^/2)    ^ 

It  follows  that 

a-\  liCjl=  0(|U^/2t"'||+a-2  lU^t'il) 

=   (|ka/2t'"ll  +  a-Vle  lUJI) 
a-\   lUJI  =   0(|U^/2*t**ll  +  hV2  lKa/2t"'ll+  a-2  lU^ijl) 

=   0(|K^/J  ||+hV2  |ka/2t'"  i  +  ea-5  IICJI  )     . 

It  follows  then,    using   (4.4)    for  sufficiently   small   e     that 


o 


(4.16)  |U^/2t'"l|  =  0(£a-V^  IKall) 

(4.17)  II  Vs'^'"  =  °^^^'^  I'^a"^  • 


This  in  turn  implies  the  polntwise  estimates 
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(^.18)       1^2^''  "  0(eh-^)  ,    Ka/2^1  =  OCea"^)  . 

Continuing  in  this  fashion  we  get  analogous  estimates  for  all 
the  higher  derivatives.   V/e  find  from  the  differential  equations 
(2.12),  (2.14)  and  (4.13a,b),  replacing  t  by  t 

IKaA^""!!  =  0(|U^//t"ll  +11C^/^P""I1  +  IUaA^"ll) 

=  0(aV2  lU^/^P  'll+ah-1  !|C^/^p""||+ah-l  I!  C^/^q"  II 
+  a-V^  iU^tll) 

llC^z/t'll    =   0(ah-l  1U^//P"|1+   llC3/4P""ll  +  IIC^/4q"ll+a-2  ||  c/t  1| )  . 

Here,  putting  M2  =  Max  HCa/^t"!!, 

II  ^a//'""-^  II  V^^"  II  =   0(IIV4^t""lUlU^/^t't'"|| 

+  IU,/,t"2|Kl|C^/^t'V|UllC3,/,f^||) 
=   0(e|lC^//'"ll+M^   IIC^/2t'"ll  +  M2lKat"|| 

=  0(e  |U^^i,t""ll+ea-V\||Cll  +  eW2  ||C^||) 
Ik^/4*P'||    =  0(|U3^/^tt""|klU^/,,  tt'"||+||C^/^t't||  +  |U^/^t"t|l) 

=   0(e  IKa/^t-'ll  +m^  IK3^/2^'"il   +M^  |U^/2  '*"  H  +^2  H^a'^'l'^ 
=  0(e  |Ua/,,t""|l  +£a-%  IKall  +  e'^a-V^  lUalD    • 

Using  also   that  by  Sobolev's  inequality   (3.13) 


*»■■    i: 


■>}  ^'  a!     -^T 


"  '.i 


'\      -U." 


"  li  .."j.\i'-  }h".twO 


*  I- 


I      r 


\\V 


;•;    ■-:•'    >  :   j": 
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MgllCjl  =  0(a2|U^/^t""l|+  ||C^t"||) 


v;e  find  that 

II  W4^"""  =  0(£a-Vl  lUJI)  ,   II  Vu't'll  =  0(ea-'^  \\zj) 
and  then  that 
(4.19)    IV^t"!  =  0(£a-Vl)  ,    Ua/4t  I  =  0(ea-2)  . 

For  the  higher  derivatives  estimates  for  the  norms  are  obtained 
similarly,  even  without  recourse  to  Sobolev's  lemma,  and  then 
pointwise  estimates  with  Sobolev's  lemma.   One  finds  generally  at 
a  point  with  x,  =  x, ,  Xg  =  Xp,  where  C  =  1  that  for  any  n-th 
derivative 

It^-^^l  =  OCea^-V^) 


with  the  sharper  estimate 


It^^^l  =  Olea-^) 


if  only  differentiations  with  respect  to  x^  and  Xp  are  involved. 
As  a  consequence  we  have  generally 
differential  equations  the  estimate 


As  a  consequence  we  have  generally  for  the  non-linear  terms  in  the 


This  enables  us  to  extend  the  better  estimate  for  jt    |  to  some  of 
the  n-th  derivatives  of  t  that  involve  x^-differentiatlons.   V/e  find 
from  the  differential  equations  (2.14)  for  h  =  a,  k  =  P  that 

It^^'l  =  0(ea-") 


.  M  «  t 


.-i^iM    u.ii'i  e»; 


'  {! .  /  i!  ""?<"'  s:r 


M     I    « 


.     >• 


S^-)o 


!>■■ 


-i(;^    rf=.rj    iHai 


(tx.n 


.^jj.. 


t  ;.'l-i.(  ,C:  ; 
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whenever  an  even  number  of  x-,-differentiations  are  involved,  and 
then  from  (2.12)  that 

It^^h  =  0(sa-") 

when  an  odd  number  of  x^-differentiations  are  involved.   Then  again 
from  (2.12) 

It^^^l  =  0(ea-^) 

for  an  even  number  of  x-,-differentiations. 

Consider  generally  an  n-th  derivative  of  a  general  component 
of  t: 

D^t,^  ,   dJ  =  dJi  D^^  D^5  ,    ^r  =  J^  ' 
where 

Denote  by  v  the  number  of  subscripts  of  t.,  that  have  the  value  5, 
so  that  V  takes  the  values  0,  1,  2.   v/e  have  then  at  the  points 
with  X  =  X  that 

(4.20  l^'^^^ik'  "  Oiea-'^-^'h^) 

where 

p,  =  0   for   j,  +  v   even  or   j^  =  0 
(4.21) 

p,  =  -1   for   J^^^   ^^^  ^^'^        J3  ^  0  • 

We  can  use  the  boundary  conditions  (4.2)  to  improve  this  result. 
Since  for  any  j,,  jg,  k 

^1^  ^2^  ^3k  =  °   ^'"'       ^3  ==  -^ 


:-'^,;«:fii 


''t       .V 


*  • »  ' 


I  .;^.  V 


wjo 
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it  follows  that 

J-i   3o 
^  ^2  ^3  bk  =  ^ 

for  some  x,  with  |x^|  <  h.   Then  for  any  x^  with  |x-,|  <   h 

^i'  ^2'  °5  bk  =  °^h  ^5^  ^^  ^2"  °|  ^5k^ 
and  also 

^1   ^2   bk  =  °<^  ^S!!  ^   V  5  bk^  • 


It  follows  that 

r 


5 


p.  =  2   for   J   =  0  ,    V  =  2 


(4«22)     '   <  M.  =  1   for   j,  =  1  ,    v  =  2 


(i  =  1   for   j,  =  0  ,    V  =  1  . 


In  particular  we  find  for  Ji  =  Jp  "^  ^  that 

(4.23)       It^^l  =  0(eh^a-^)  ,    |t.^|  =  Oleha"^)  . 


5.  Estimates  for  the  displacement  gradients 

Relations  (4.21),  (4.22)  summarize  our  results  for  the  orders 
of  magnitudes  of  the  derivatives  of  the  various  components  of  the 
stress  tensor  t.   They  hold  under  the  assumption  that  there  exists 
a  cylinder  of  radius  a  about  the  point  in  question  in  which  the 
differential  equations,  the  boundary  conditions,  and  the  estimate 
(4.3)  for  t  hold.   Moreover  e,  h,  a  shall  satisfy  the  inequalities 
(4.4),  (4.4a)  where  e  is  a  certain  universal  positive  constant. 
Here  the  relation  e  <  e  is  essential  for  the  arguments  used;  it 
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has  the  effect  of  making  the  coefficients  of  the  highest  derivatives 
in  the  differential  equations  approximately  constant.   Assumption 
(4.4)  on  the  other  hand  can  be  weakened  slightly.  We  get  the  same 
results  if,  say, 

(5.1)  el4'   111.   ^'-o- 

assuming  that  the  universal  constant  e  was  less  than  1.   Indeed, 
if  e,  h,  a  satisfy  (5.1)  and  here  even  e  <   e   h^a"^  then  (4.4)  is 

satisfied  and  our  earlier  results  hold.   If  on  the  other  hand 

2  _2         *     1/2         *  '?  "> 

e  >  e^h  a   we  put  a  =  ae^'  .  Then  a  ^  a,  s  j^  e^h  (a*)"  ,  and 

h(a*)"  =  ha'-'-e"-'-/^  <  ^''^^'^  ^   1.   Our  previous  results  (4.21), 

(4.22)  then  apply  with  a  replaced  by  a*.   Since  however  a  and  a* 

only  differ  by  a  universal  factor  we  can  again  write  a  for  a*, 

absorbing  the  factors  into  the  symbol  "0". 

The  displacement  gradients  p  ,  =  x.  ,  are  connected  with  the 

g.,  by  relations  (1.4).  Here  an  overall  restriction  on  the  length 

of  the  plate  considered  becomes  necessary  to  prevent  rotations  from 

varying  too  wildly.  To  simplify  matters  our  "plate"  in  the 

unstrained  state  shall  be  a  circular  cylinder 

x^+Xg  <  4a^  ,    |x^!  ^  h 

in  v/hich  the  differential  equations  (2.12),  (2.14),  the  boundary 
conditions  (4.2)  and  the  estimate  lt|  <  £  shall  hold.  Moreover  e, 
a,  h  shall  be  restricted  by  the  inequalities  (5«1)«   In  the  smaller 
cylinder 

(5.2)  x^  +X2  1  a^  ,•    Ix^l  <  h 


■'u' 


^    ■u 
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the  stress  derivatives  will  then  satisfy  the  inequalities  (4,21), 
(4,22).   For  the  following  estimates  of  the  displacement  deriva- 
tives we  shall  restrict  ourselves  to  the  smaller  cylinder  (5.2). 

We  first  of  all  apply  a  rigid  motion  to  the  deformed  plate  to 
remove  the  rotation  at  the  origin.  Then  the  matrix  p(0)  =  (p.,(0)) 
is  symmetric.   Since  the  square  of  this  matrix  is  the  matrix  with 
components  g.,  =  ^-iu;+ ^^iv  =  5.,  +  0(e)  we  have  then 

(5.3)  '^  9r        Pik^O^  =  5.j^-i-0(e). 
By  (1,6)  for  any  x  in  (5-2) 

PjiPji  =  Sii  =  0(1)  ; 

moreover  the  square  of  the  determinant  of  the  p..  equals  the 
determinant  of  the  g.,  which  has  a  value  l+0(e).  Hence  also  the 
matrix  reciprocal  to  that  of  the  p^,  is  bounded: 

p^^  =  0(1)  ,   pi^  =  0(1)  . 
It  follows  from  (1.12)  that 

(5.4)  p'  =  O(g')  =  0(£h"^)  , 


(5.5)     Pii,(x)  =  p.j^(0)  +  0((a+h)|p'|)  =  5^^+0(eah-^) 


Then 

By  (1.6),  (5.5),  (5.1) 

(5.6)     (Pik-  ^Ik^  -^(Pki-  W  =  ^Sik-^ik^  -  (Pji-  ^ji^^Pjk  -  V^ 

=  0(e+  e^a^h'^)  =  0(e)  , 
In  particular 


::,  .  &n     &\ 


^d3  r:zSz  ', -^aH 
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(5.7)  Pll.P22'P53  "  ^"^^^^^  • 
By  (1.4)  then 

(5.8)  P^^-^ik  =  -(Pik-^k^  +  O^^Vh-^)  =  0(eah-^)  , 

(5.9)  p-'-^,P^^,P^^  =  l+0(e)  . 
By  (1.12),  (5.8) 

=  p^l[P,  2  a  ]+p^^[3,  2  a  ] 
=  0(|t|  +  eah"^lt'| )  =  0(£a"^) 
while,  of  course,  p-j^g  ,  =  0(eh"  )  by  (5.4).  Then 

^12  "  Pia^'^^  +  ^^^IPlo' "^^1^12  5'^  "  °^^^  • 
Combining  with  (5.7)  we  have  then 

(5.10)  Pap-^'   p^^-l  =  0(e).       .,=-?=.s^e 

The  first  of  the  v.Karman  equations  is  based  on  the 
compatibility  equation  (2.l6a).   By  (2.12) 

^12,12  "  "''^22,22"*  '^25,23"  ^2k,k2 


and  also 


Using  also  that 


^12,12  ~  ""^11,11"  ^13,13"  ^Ik, Ik 


^,11  "^^,22  -  ^aa,pp'^^33,PP  ' 


relation  (2.l6a)  can  be  written  in  the  form 
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t5-ll)  *^aa,Pp-^Paa,pp+(l  +  ^)ta5,a3 


=    {l+cr){q  _   p  )  +   a-t  +P 

Here  by  (4.22) 

(5.11a)  *33,pp  =  0(£h^a"^)    . 

Similarly  by  (2.13) 

(5.11b)   PaP,75  "  0(|t  V|+|t^|)  =  0(e2a-^)  =  0(eh^a"^)  . 

Using  (4.22)  we  also  find  that 

P_  ^=(t  t,-  -^^  t.  .t  ,)  .   =^e^a"^)  =  0(£h^a"^)  . 
a3,a3     ar  r3   1 +cr  ii  a3  ,a3 

Finally  by  (2.15) 

^1212  =  2[a,  1  2]{^y^  -2[a,  1  Hj/g^ 

=  2[a,  1  2]s^^[b,  1  2]-  2[a,  1  l]s^^[b,  2  2]  . 

—  1       ab 
Since  each  3-index  symbol  is  of  order  0(eh~  )  and  g   =  ^ab"*"^^^^ 

we  have 

^1212  "  2[a,12][a,l  2]  -  2[a,l  l][a,2  2]  +  0(e^a"^) 

=  2([3,12][3,12]  -  [3,ll][3,22])  +  0(e2a-^)  . 
Hence  (5.11)  yields 

(5.12)   t^ct,PP  +  Paa,Pp-^(l+-^)^a3,a3 

=  2(1+  T)([3,12][3,12]  -  [3,ll][3,22])  +  0(£h^a-')  , 

We  now  introduce  the  quantity  w  =  x^(x^,X2,0)  which  gives  the 
normal  displacement  of  the  "middle  surface"  x^  -  0.   We  have 


^+5 


^^aP  =  P3a,P^^1^^2'0) 


Now  by   (1.12),    (5.8) 


P3a,p  "  P^^[k,a3]   =   [  3,ctP]  +  0(  e^h"^)    . 
Since  also  p,^  ^  =  0(eh'-'-)   by   (5.^),   we  find  that 

[3,12][5,12]  -  [5,11][3,22]    =   P5i^2P31, 2-^31, 1^32, 2"^°^^^^"^^    ' 
Moreover,   by   ('4.21) 


[5,aP]^3  =  i  (S3^,P3''^3P,a3-SaP,33^ 


1  ^ 

=  2   ^^3a,P3  "^^3^,03  " ^ap , 33  " TTr ^oS  *il,33^ 

=  0(ea"^) 


and  hence 

lit.  Toir-z  Toi      r^  Tiir^  oo'w 

It  follows  that 


([3,12][3,12]  -[3,11][3,22])    ,  =  0(£W^)    . 


[3,12][3,12]  -[3,11][3,22] 


=    ([3,12][3,12]  -[3,11][3,22])^  ^Q+0(e2a-2) 

3 

=  w^^2^l2-^ll^22-^°(^^^"^)  • 


Thus   (5.12)   can  be  written  as 


2   -4 
=   2(1  +t)(w  ^gW  ^2  -^  11^  22^  +0(£h  a      ) 
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We  introduce  nov;  the  average  stresses  In  the  tangential 
directions  over  the  thickness: 

-h 

As  a  consequence  of  the  differential  equations  (2.12)  and  of  the 

boundary  conditions  t  ,  =  0,  which  also  imply  that  P  ,  =  0,  we 

obtain 

h 

-h 

The  "to  can  then  be  expressed  with  the  help  of  an  "Airy  function" 
0  =  0(x-,,X2)  in  the  form 

Hi  "  ^,22  '    "^.12  "  "^,12  '        "^.22  "  ^,11  • 
Integrating  (5»15)  over  the  thickness  and  using  that  t  -,   =  0  on 

v-X  ^  ^  tx 

the  boundary,  we  obtain  immediately  the  first  equation  of  v.  Karman 
in  the  form 

^5-15)     0,^«pp  =  2(l+o')(w^3_2^l2-^ll^22)  +  o(^^^^"'')  • 

In  order  to  obtain  the  second  equation  we  have  to  calculate 
the  fourth  derivatives  of  w.  By  (1.12),  (1.13) 


k3 
^3a,a  -  ^ 

and  then 


P3a,a  =  P'^tk,aa] 


5\1 
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P3a,app  =  P^W^"[i,sP][j,rP][k,aa] 

=  -[k,5P]^p[k,cia]  -  [k,3P][k,aa]   p+[5.cta]   ^^ +0(  e^a'^h"^) 
"  "  2  TT^  ^aa,3PP  "^  T^  ^33,3PP+^(^ha'   ) 


Since 


irif  ^au.jpp-^"'*^"*)  • 


*aa,3a35  =  °<^^'*' 


we  can  replace  t   ^gg  by  its  value  at  x^  =  0  without  changing  the 
order  of  the  error  term.   It  follows  that  for  any  x. 


-2  ^aa^^=TT^taa,3^^■*•°^'^"''^^  ' 


Then 


h 

-h 

h 


-h 
h 
i-^     I     X.,  t  ra   dx-,  +0(eh^a~^^)  . 
-h 


fo 


41  *~  i'»*sJx-^.Ji 


t-^-xI 


v.—r   ^  -.':     *■ 


I      n 


('■  v"  J:.;. 
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By  (5.13),  (5.11b)  we  have  then 

h 


(5.16)   |h5w^^^pp=  (1-^)/  ^3*a3,a5^^3-*-°^'^''^"'') 


-h 


h 
-h 


Here 


h 
-h 


=  (1-qr)  /  P^j^  j^  dx  +  0(ehV^)  . 


^3k,k  "  ^^5r^rk"r+T  ^ii^3k^k"T  ^^ir^ir'TTo^  ^ii^rr^3 
'  "^iS^rk,k  +TTo^  ^ii^3k,k  ^  ^3r,k^rk'  TT?  ^ll,k^3k 

■  2  ^ir^ir,3  2  1  +  Q'^li^rr,3 
=  -I  WaP,34r?^^aa^3P,3-'^^^^^^"^^ 
=  -  I  tapSap,3-^°(^^^^"^^  =  [3,aP]t^p+0(e2ha-2) 

Substitution  into  (5,16)  yields  then  the  second  equation  of 
v.Karman 

(5-17)   w^^^pp  =  |(l-^)h-2(w^^^0^22-2w^^2^^^2+^22^11^+°(^^^'''^ 

In  the  customary  notation  one  normalizes  the  Airy  function 
differently,  introducing  as  Airy  function 

X=   2(1+  r)  ^ 


>  if -I 
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so  that  X   has  the  same  relation  to  the  ordinary  stresses  as  our  0 
has  to  the  dimensionless  stresses  t.,  .  The  two  equations  (5»15)* 
(5.17)  then  take  the  form 

(5.18a)        \  [? ]C   =  w  ^gW  3_2  "^  11"  22  +  °^^^^^'  ^^  * 

(5.18b)   A^w  ^\   (w^llX,22-2^12/^,12+^22/,ll^■'°^'^^''^^ 
where  N  =  Ehv3(l  -  cr  )  is  the  bending  stiffness  and  where 

A  =  — ^  +  — ^  . 

Sx^    Sx2 
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